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Abstract 

We study scattering of waves by black holes. Solving a massless scalar field with a point source 
in the Schwarzschild spacetime, waves scattered by the black hole is obtained numerically. We then 
ff^ reconstruct images of the black hole from scattered wave data for specified scattering angles. For 

^ the forward and the backward directions, obtained wave optical images of black holes show rings 

^ that correspond to the black hole glories associated with existence of the unstable circular photon 

^ orbit in the Schwarzschild spacetime. 

^ PACS numbers: 04.20.-q, 42.25.Fx 
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I. INTRODUCTION 



Wave scattering and resulting diffraction effect are related to familiar phenomena such as 
rainbows and glories. These phenomena are caused by scattering of light by water droplet 
in atmosphere and analysis based on the wave optics is required to understand its formation 
and quantitative nature [Ij. The methodology and the formalism of the wave scattering 
problem are applied to wide fields in physics for the purpose of obtaining the information 
of the scatterer by analyzing the date of the scattering wave. For the gravitational physics, 
the wave related properties of black hole spacetimes have been investigated by many re- 
searchers for more than 40 years to obtain physical nature of black holes [2]-[9]. T. Regge 
and J. A. Wheeler [2j investigated the wave equation for gravitational perturbations in 
black hole spacetimes and have shown that the Schwarzschild black hole is stable against 
the perturbations. For the spinning Kerr black hole, it was shown that the incident wave 
is amplified if the wave satisfies an appropriate condition. This phenomena is called super 
radiance and peculiar to the Kerr spacetime accompanying with the dragging effect [lOj . 
The quasi-normal oscillation of black holes is also wave related property and this normal 
mode is characterized by parameters of black holes fTO]. As a straightforward application 
of the wave scattering theory, cross sections of black holes are obtained from the scattering 
waves and the nature of the black hole geometry is discussed from the view point of waves. 
Especially, peculiar to the wave scattering by the black hole, the diffraction effect for the 
backward scattered wave becomes significant and this leads to the phenomena so called the 
black hole glories ^ [5] . The main purpose of the wave scattering problem by black holes 
is theoretical understanding of the physics of black holes and wave propagation in curved 
spacetimes. For astrophysical black holes, it is unlikely that wave effects such as diffraction 
can be observed using current technology, but there is a possibility that future technology 
enable us to study interference effects in gravitationally lensed waves. 

The treatment of the conventional wave scattering by black holes mainly concerns analysis 
of the scattering amplitude and the cross section. In this paper, we consider wave scattering 
by black holes from the view point of image formation. Our main motivation is to investigate 
''images" of black holes formed by the incident wave using wave optics [HI [12]. In the 
geometric optics, "images" of black holes are obtained by solving null geodesies (light rays) 
in black hole spacetimes and we can draw images of a black hole as "black hole shadows" ; 
if we assume a light source behind a black hole is larger than the angular size of the black 
hole, a distant observer see a dark spot that corresponds to the apparent image of the 
black hole. We expect to obtain the shadow image of the black hole in the framework of 
wave scattering problems. However, obtaining images using waves is not trivial task at first 
sight; as is known, the scattering amplitude at the observer shows diffraction pattern due to 
interference between scattered waves. However, this diffraction pattern is not itself images 
of the black hole. The theory of image formation in wave optics tells that we must decode 
scattered waves using imaging device to reconstruct images [13j. The scattered waves by a 
black hole contains information of the black hole geometry and it is possible to reconstruct 
images of the black hole. For this purpose, we introduce a convex lens to the configuration of 
the standard wave scattering by black holes and aim to reconstruct images from scattering 
waves. 

In this paper, we consider scattering of massless scalar waves by a Schwarzschild black 
hole. The massless scalar waves are adopted as the benchmark treatment for wave scattering 
problems by black holes and we do not consider polarization degrees of freedom that is 
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necessary for the electro magnetic waves and gravitational waves. We investigate images of 
the black hole using wave optics. As the source of incident waves, we prepare a point source 
of the wave placed near the black hole. This models the astrophysical black holes with 
active emission regions such as accretion disks around them. We solve the wave equation 
for the massless scalar field (Helmholtz equation) using the finite difference method. Then, 
we reconstruct images from scattering wave data and investigate characteristics of the black 
hole spacetime appearing in images. This paper is organized as follows. In Sec. II, we 
introduce basic equations and our numerical method for solving the scalar wave equation. 
In Sec. Ill, we explain the method of images reconstruction from scattering waves. In Sec. IV, 
we present our numerical results and Sec. V is devoted to summary. We use units in which 
c — h — G — 1 throughout this paper. 



II. WAVE SCATTERING IN BLACK HOLE SPACETIMES 



We solve the wave equation for the massless scalar field in the Schwarzshild spacetime 
numerically. For this purpose, we rewrite the Klein-Gordon equation for the scalar field 
to the Helmholtz equation assuming that the wave field is monochromatic. Then by finite 
differentiating the equation, we obtain the numerical solution of the wave scattering problem. 



A. Scalar wave equation in Schwarzschild spacetime 

For the massless scalar field $, the scalar wave equation in a curved spacetime is 

1 
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where is a source term of the wave. For the Schwarzschild spacetime, the metric is 
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and the wave equation in the Schwarzschild spacetime is 
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We consider the solution of this equation corresponding to the wave scattering problem by 
the black hole; the wave is emitted by a source placed at r = rs, ^ = — tt (on the z axis) and 
scattered by the black hole. The scattered wave is observed by a detector placed at r = robs 
(Fig. 1). We do not consider the motion of the source in this paper. The location of the 
observer is not restricted on the z axis. We assume that the wave is monochromatic with the 
angular frequency lj. This angular frequency is defined with respect to the Schwarzschild 
coordinate time t. 
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detector's coordinates 



FIG. 1: Configuration of the wave scattering problem by a black hole. We introduce spatial 
coordinates (x^y^z) = (r sin 9 sin (j)^ r sin 9 cos (j)^r cos 9). A point source is placed at (x^y^z) = 
(0, 0, — rs). The scattered waves are observed by a detector at (x, y^ z) — (robs sin^o? 0, robs cos^o)? 
where corresponds to the scattering angle. 

Because the Schwarzschild geometry is spherically symmetric, the source and the black 
hole system of our scattering problem has an axial symmetry. Thus we can impose the 
axial symmetry for the wave configuration: 5^^$ = 0, d^j^S = 0. As we assume that the 
time dependence of the scalar wave is $ oc e~^^^, the wave equation can be written as the 
following Helmholtz equation: 
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where $ = r$ and the tortoise coordinate is introduced by x 
denotes the two dimensional delta function. 



2Mln(r/2M- 1). 5^ 



B. boundary conditions 

We must impose boundary conditions to obtain solutions of our wave equation. We first 
consider the boundary conditions at r = 2M and r = oc. By the definition of black holes, 
the wave is purely ingoing at the horizon r = 2M. Near the horizon, the Helmholtz equation 
becomes 

dl^ + uj^^ = 0. (5) 

This equation contains the derivative with respect to the radial coordinate only, and the wave 
propagates perpendicular to the horizon. Taking into account the purely ingoing condition 
for the wave at the horizon, the wave goes through the horizon perpendicular to it. By 
restoring the time dependence, this boundary condition can be expressed as 



-oo 



= 0. (6) 



On the other hand, the wave is purely outgoing at the spatial infinity because the source of 
wave is located at the finite distance from the black hole and no wave is coming from outside 
of the numerical box. Thus we assume the wave behaves approximately as a spherical wave 
at the outer boundary. This boundary condition can be expressed as 

dj + dt^) _ =0. (7) 
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Next, we consider the boundary condition associated with the symmetry of the system. 
According to the axial symmetry of the scalar field configuration, $ should be an even 
function with respect to the angle variable 0: $(^) = $(— ^),0 < ^ < vr. Hence, the 
boundary condition for the wave on z axis is 



= 0. 



(8) 



In Eq. ([4|), the first order derivative term cot 9do^ is contained and its coefficient diverges 
at = 0,7r. Around = 0, considering the boundary condition (Is]), 
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9 

where ' represents the derivative with respect to 9. Hence, 

cot9^'{9) ^^"{0) + O{9) 



(10) 



and this term is regular at ^ = 0. In the similar manner, we can see that the equation ^ 
is regular at ^ = tt. After all, at ^ = 0, tt, Eq. ^ reduces to 
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We adopt this equation as the equation to obtain the wave field on the z axis {6 = 0, tt). 



C. finite difference method 

To solve the Helmholtz equation numerically, we should rewrite the equation into the 
corresponding difference equation. For our numerical calculation, we divide the considering 
spatial region in (x, 9) space into (A^ + 1) x (A^ + 1) homogeneous grids. The coordinate 
variables are 

X, = xo + i-A, ej=j'Ae, =0,--- , A, (12) 

where A and Aq are the spacing of the grids in x and 9 coordinates, respectively. The 
coordinates range of the numerical box is 



xo<x<xo + AA, 0<9 < NAo = tt. 



(13) 



Accordingly, the field l> at the point {xi^9j) is written as l>^j = l>(x^,0j). Then, Eq. Q in 
the finite difference form is given by 
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and Eq. (11) is rewritten as 









+ 




(-91 









hj = (15) 

where we have introduced S = Aq/ A and is defined as = r(x^). In these equations, the 
radial coordinate is measured in the unit of the black hole mass M. 

The boundary conditions for the r direction of the Helmholtz equation must be imposed 
using Eq. ([6) and Eq. Q to determine ^-ij and l>7v+i,j. These equations are conditions for 
temporal behavior of waves. Thus we have to translate these conditions for the wave field to 
that does not contain the time dependence. We first consider Eq. ([g]). This equation relates 
the wave field at the inner spatial boundary xq and As the solution of this equation is 
x) = $(t + x), the field value $(t + A^, is expressed as 

d(t + A,, = x_i + A,). (16) 

As the time dependence of the wave is assumed to be oc e~^^^^, this boundary condition 
^ provides the following relation between the field values at Xq and 

g.M.A4_^^^. ^ + ^^^^^^ _ (17) 

Thus, can be determined by 

= i_^^re-.M.^.A ^o,, (18) 

where 6t = A^/A. The boundary condition at the outer spatial boundary can be obtained 
in the similar way and we have 

^^+1'^- = i_^^_e-^Mc...A ^^.r (19) 

In our numerical analysis, we set St = I and we adopt the following equations to impose the 
boundary conditions at the inner and the outer spatial boundaries: 

= e^^'^^doj, ^N+i,j = e'^^^^Nj. (20) 
The boundary conditions on the z axis are determined by Eq. ([s]) and 

^^,-1 = ^^,N+l = ^i,N-l- (21) 

In the neighborhood of the point source, $ behaves as a spherical wave in fiat spacetime 
and we can impose the boundary condition around the point source using this behavior. We 
assume that the wave around the point source placed at xs is distributed as 

l> = A .^ ^ ^ iuj\x-xs\^ - ^ (xsin6^cos(/), xsin6^sin0, xcos6^), (22) 
\x — Xs\ 
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where a constant A represents the strength of the source. We use this relation to determine 
the wave field around the point source. In the numerical grid space, the location of the 



source is (x^g, 0). We consider the $ at points (x^g-i^ 0)' (^^s+i^ 0) (^^s' ^o)- Using (22), 

(23) 



(f) — i ^^^^ 



The distance between the point (x^g, Aq) and the source point is r^gA^i. Hence, $ at (x^g, Ao) 
is 



These relations provide the boundary condition around the point source. 



D. Solving finite diflFerence equations 

These diflFerence equations and the boundary conditions constitute the simultaneous equa- 
tions. The simultaneous equations are written in a matrix form as 



A' X = h 



(25) 



where x is the vector to be solved and the matrix A and the vector b are determined by 
the diflFerence equation. As demonstrations, let us first consider the 3x3 grids case. If we 
assume that there is no source. A, x and b are 
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where and represent the real and the imaginary part of l>, respectively. Ci, C2, Di 
and Hi are 6x6 matrices given by 
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tti = rf5^sm{MuA), 

a = S^A'^ {M^oo'^rf - 2{ri - 2)) - 2rfS^ - Anin - 2) + nS^ cos{MujA), 
di = S^A^ {M^uj^f - 2(n - 2)) - 2rfS^ - Aniri - 2), 
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ki = 2(52A2 (mW4 - 2{ri - 2)) - 44S^ - inin - 2), 



Di = diag[4ro(ro - 2),4ro(ro - 2),4ri(ri - 2),4ri(ri - 2),4r2(r2 - 2),4r2(r2 - 2)] , (29) 
Hj = diag [ro(ro - 2)(2 - cot ^j), ro(ro - 2)(2 - 5A cot ^j), ri(ri - 2)(2 - 5A cot ^j), 



n{ri - 2)(2 - 5A cot ^j), r2(r2 - 2)(2 - 5 A cot 9j),r2{r2 - 2)(2 - 5A cot 6'j) 



In the case of 4 x 4 grids, the structure of x and A becomes as follows: 
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Di diag 



4ro(ro - 2),4ro(ro - 2),4ri(ri - 2),4ri(ri - 2),4r2(r2 - 2),4r2(r2 - 2) (34) 



,4r3(r3-2),4r3(r3-2) 

i?^- = diag ro(ro - 2)(2 - 5A cot ro(ro - 2)(2 - cot ri(ri - 2)(2 - 5Acot%), 

ri(ri - 2) (2 - cot ^^^^ r2(r2 - 2) (2 - 5Acot9j),r2{r2 - 2) (2 - 5Acot%), (35) 
r3(r3 - 2) (2 - 5A cot ^^^^ r3(r3 - 2) (2 - 5 A cot%)]. 

If the number of the grids becomes larger, it is possible to predict the structure of A from 
these examples. Of course, we should modify the forms of A and b to introduce a source of 
waves along the method we have described in Sec. IIC. 

In our calculation, the region of the numerical box is —2 < x < 23 (2.03 < r < 20.5), < 
9 < IT and we divide this region into 1001 x 1001 grids. We used the Mathematica to obtain 



the solution of the linear system (25). 



III. IMAGE FORMATION IN WAVE OPTICS 



Most of previous works on wave scattering by black holes aim to obtain the scattering 
amplitude or the differential cross section to investigate wave effects. On the other hand, 
we want to consider images constructed from the scattered waves by black holes. By using 
images, it is also possible to discuss the semi-classical nature of the wave scattering by black 
holes. We have introduced the method to reconstruct images from waves in the gravitational 
lensing system [HI [12] and obtained images of gravitational lensing in wave optics. We 
apply the same method to the wave scattering by black holes. In the wave optics, the 
image formation can be understood as diffraction effect of waves. For this purpose, image 
formation devices such as a convex lens are introduced. 
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FIG. 2: An image formation system with a convex lens. Incident waves are transformed by the 
lens and an image of a source object appears on the image plane. This image corresponds to the 
Fourier transformation of the incident wave. 



In our analysis, we adopt a thin convex lens. Let us $(X) be the incident wave in front 
of the lens and ^^{X) is the transmitted wave by the lens. X = (X, Y) denotes the location 
of a point on the lens plane Z = 0. These two wave functions are connected by the following 
relation 

$t(X) = e-^"^$(X), (36) 
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where / represents the focal length of the lens. For a point source placed at Z = — / (front 
focal plane), the wave in front of the lens is 



where we assumed |X| <C /. Hence the transmitted wave becomes 

$t(^) = e'""^ 

and its phase becomes independent of (X^Y). This means that the transmitted wave is 
the plane wave. Thus the convex lens converts a spherical wave front to a plane wave 
front. Now we consider the wave amplitude on the focal plane Z = f (image plane). Using 
the Fresnel-Kirchoff diffraction formula [13j, the wave on the image plane is given by the 
following diffraction integral: 



$i(Xi) oc / rf2X$(X)e-^"^ X 

'\X\<d 



(37) 



where £ is path length between a point on the lens plane and a point on the image plane 
(see Fig. 2) and d is a radius of the lens. Using I = J\X -Xi\^ + P ^ f + \X - Xi\^/2f, 



$i(Xi) oc / dX <!?{X) e-f^^^-^\ 

'\X\<d 



(38) 



Hence $i is the Fourier transformation of the incident wave $. 

Let us consider the situation that the incident wave can be written as the WKB form 



$(X) = yle^'"^(^) =^exp 



iio (^S{X,) + ^S"{X,){X - X,f + 



(39) 



where = (X*, Y^) denotes the location of the point that the null geodesies intersects with 
Z = plane. This point corresponds to the classical path of the null ray obtained as the 
saddle point of the action S\X^) = 0. The classical path is a null geodesies connecting 
a source and a point on the lens plane. For this form of waves, the wave amplitude on the 



image plane (38) becomes 
$l(Xi) cx 



(fx exp 



\X\<d 
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S"X,+Xi/f 





(40) 



We have assumed \(juS^^(f\ <C 27r and omitted the quadratic term |Xp in the exponent to 
evaluate the integral at the last step of the calculation; this is the condition for the Fraunhofer 
diffraction. On the dimensional analysis, S" ^-^ 1/r where r is the distance between the black 
hole and the observer. As <C r, this condition is easily satisfied. For ujd^ 1 (geometric 
optics limit), we obtain 

$i(Xi)oc52fXi + /5"X,). (41) 



10 



Thus we recover a point image of the point source on the image plane; the location of the 
image corresponds to the location determined by the null geodesies. If we does not take 
geometric optics limit, the image acquires diffraction effects and the image of the point 
source has a fine size called the Airy disk. We apply this model of image formation to 
scattering problems by black holes and obtain images of black holes (Fig. 1). For a point 
(X, Y) on the lens plane Z = 0, the angle between this point and the z axis is given by 



cos 9 = — sin^o + A/ l— ( — 
r \ \ r 



Y 



cos 9f) . 



(42) 



We use this relation to evaluate the Fourier integral (38) to obtain images. 



As an example of image reconstruction from scattering waves, we present the gravitational 
lensing by a point source. Assuming that the gravitational field is weak, the wave equation 
for the massless scalar field reduces to the Newtonian form [4J 



4Ma;^ 



$N = 0, 



(43) 



where is the flat space Laplacian and the problem is equivalent to that of the Coulomb 
scattering in quantum mechanics. The exact wave function for the scattering problem with 
a plane wave incident from infinity is given by 



$j^(^o) = e'^'^'^T(l - 2icoM)e" 



iFi{2iuM, l;iur{l - cos 6*0)). 



(44) 



The wave amplitude at r = 20M is shown in Fig. 3. For small scattering angle, oscil- 
latory behavior due to diffraction between the incident wave and the scattered wave can 
be observed. The behavior of the scattering amplitude for small scattering angle (forward 
direction) can be obtained by the asymptotic behavior of the confluent geometric function: 



^^{9o) (X Jo{2ujVm^9o) for l^ol 



< 1. 



(45) 



The Fourier transformation (38) for this function results in the function with a peak at the 



angular radius 9^ = \/4:M/r^ that corresponds to the Einstein ring. On the other hand, for 
large value of 9q^ the diffraction effect becomes small and a ring image does not appear in 
the backward direction. 
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FIG. 3: The wave amplitude |$n| for Muj = 12 at r = 20M. 
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The reconstructed images from the scattering wave using (38) are shown in Fig. 4. For 
^0 = 0, we obtain an image of the Einstein ring with angular radius 6e- For 7^ 0, double 
images of the point source appear. These results reproduce images of the gravitational 
lensing obtained by the ray tracing method using the lens equation. 




FIG. 4: The left panel is the image for = and this is the Einstein ring. The angular radius of 
the ring is given by ^e- The right panel is the image for = ^/4 and double images of the source 
appear. These images are obtained for ujM —12^r — 20M, rf/r = 0.2. 



12 



IV. NUMERICAL RESULTS OF WAVE SCATTERING BY BLACK HOLE 



We numerically solved the massless field in the Schwarzschild spacetime for the angular 
frequencies Muj = 2, 12, 24 and the source positions rs = 2.5M, 3M, 6M, 15M. The detector 
is located at robs = 20M. The numerical grid size is 1001 x 1001. We have also done the 
calculation with the grid size 501 x 501 to check the validity of our numerical results. In 
this paper, we only present the result for Muj = 2, 12 and rs = 2.5M, 6M. 

A. Muj 12 case 

We first present the result for Muj = 12 case. The source position is rs = 6M and 
moderately far from the black hole. This setting corresponds to the standard analysis of 
wave scattering by black holes that the source is placed at spatial infinity and the incident 
wave is treated as a plane wave. At the location of the observer, we do not distinguish 
the incident waves and the scattered waves. Hence the detector receives both the incident 
waves and the scattering waves, and their superposed wave amplitude is obtained. As the 
result, despite of the long range Coulomb like nature of the gravitational field, "scattering 
amplitude" |$obs| remains finite even for small scattering angles (see Fig. 6). The situation 
is the same as the gravitational lensing presented in the last section. The standard handling 
of the scattering problem extracts purely scattered wave and results in diverging scattering 
amplitude for small scattering angle (forward direction) due to the long range nature of the 
gravitational force. 




-4 -2 2 4 -4 -2 2 4 



FIG. 5: Spatial distribution of $ for Muj = 12 on {z,x) plane. Left panel: Re[$]. Right panel: 
Im[$]. The point source is placed at (z, x) = (— 6M, 0). 

Fig. 5 shows the distribution of $ on (z, x) plane. We can observe that the wavelength 
of the incident wave becomes shorter near the black hole horizon due to gravitational blue 
shift. Along the z axis (z > 0), the amplitude of the wave is enhanced due to constructive 
interference of scattered waves. It is also possible to observe the characteristic wave pattern 
formed around the unstable orbit r = 3M. 
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FIG. 6: The observed wave amplitude at robs = 20M. Muj = 12, rs = 6M. 



As already mentioned, the wave amplitude remains finite for 0. The wave amplitude 
has oscillatory behavior for all values of the scattering angle. This is contrasted with the 
scattering amplitude for the gravitational lensing. Except for 0, the wave amplitude 
increases as the scattering angle increases. As the distance from the point source to the 
observer varies depending on the scattering angle, the observed wave amplitude increases for 
the backward direction because the distance between the source and the observer decreases. 
The oscillatory behavior of the wave amplitude at the forward and the backward directions 
corresponds to so called glory effect [3l [5] and it is possible to identify the glory in images 
in our analysis. For black holes, glories arise because a wave can be defiected through an 
angle greater than tt and the glory scattering is associated with the unstable photon orbit 
at r = 3M. 

The reconstructed images from the scattering wave are shown in Fig. 7. We have applied 
the formula ([38]) with d/vohs = 0.5 to obtain these images. For Oq = 0, the image is a ring 
that corresponds to the forward glory. As is known, the scattering amplitude for the glory 
is given by ^ ^ for the forward and the backward direction 



(46) 



where £g = 3\/3 Mu corresponds to the critical impact parameter that the incident nuU ray 



can escape to infinity. In terms of images obtained by the Fourier transformation of (46), 
the apparent angular radius of the ring is given by 



3^/3M 



(47) 



For the backward direction = ^r, we cannot identify the ring image well because it is 
drown in the wave of the point source. But in the plot of the intensity distribution (Fig. 8), 
we can identify small dips of the intensity at the radius corresponding to the ring of the 
backward glory. 
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FIG. 7: Images of black holes reconstructed from scattering waves {Muj — 12, rs = 6M). From the 
top left panel to the bottom right panel, the scattering angles are — 0, 7r/4, 7r/2, 37r/4, tt. 




FIG. 8: The intensity distribution of images. Left panel: = (forward). Right panel: — tt 
(backward) . 

In the geometric optics limit, images of black holes can be obtained by solving null 
geodesies. For the observer at = 0, the primary null rays, which are deflected by the 
black hole but do not go around it, result in the Einstein ring. The secondary and the 
higher degrees of null rays that go around the black hole many times also form ring images 
with smaller angular radius compared to the Einstein ring. The left panel in Fig. 9 shows an 
example of the primary and the secondary null rays connecting the source and the observer 
at ^0 = 0. The right panel in Fig. 9 shows the apparent angular radius of the Einstein ring 
and the ring by the secondary rays as the function of the source position rs- For rs = 6M, 
the ratio of angular sizes of these two rings are ^1.2. 
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FIG. 9: Left panel: the primary and the secondary nuh rays connecting the source at rs = 6M 
and the observer at robs — 20M. These rays constitute ring images. Right panel: The apparent 
angular radius of ring images for different location of the source. The blue line corresponds to the 
primary rays (Einstein ring) and the red line corresponds to the secondary rays. 

We expect to observe these double ring structure in our image reconstructed from the scat- 
tering waves. However, as are shown in Fig. 7 and Fig. 8, it is not possible to identify this 
structure because it is buried in the diffraction pattern appeared in the image. By using 
higher frequency waves, identification of the double ring structure of the images will be suc- 
ceed. For this purpose, we present the image using the wave with Muj = 24 (Fig. 10). The 
highest peak of the intensity corresponds to the unstable orbit and it is possible to observe a 
broad second peak outside of it. The measured radius of the second peak is about 1.2 times 
larger than that of the unstable orbit and we can conclude that the second peak corresponds 
to the Einstein ring. 




FIG. 10: Left panel: image for = using the wave with Mcj = 24, d/robs = 0.6. Right panel: 
intensity distribution of the image. 

We then consider the case that the source is near the black hole rs = 2.5M; in this case, 
the source is placed inside the unstable orbit 3M. We do not observe significant difference 
of the scattering behavior of waves compared to rs = 6M case. 
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FIG. 11: Spatial distribution of $ on {z,x) plane for Mcj = 12. Left panel: Re[$]. Right panel: 
Im[$]. The point source is located at {z,x) = (— 2.5M, 0). 
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FIG. 12: The scattering amplitude at robs = 20M. Muj = 12, rs = 2.5M. 

In this case, it is possible to observe a ring for the backward direction Oq = n in image 
(Fig. 13) because the distance between the source and the observer is larger than rs = 6M 
case and the amplitude of the wave that directly reaches the observer reduces. This enables 
us to observe the backward glory as a ring image. 
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FIG. 13: Images of black holes reconstructed from scattering waves. From the top left panel to the 
bottom right panel, the scattering angles are — 0, 7r/4, 7r/2, 37r/4, tt. 
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FIG. 14: The intensity distribution of images. Left panel: 0q = 0. Right panel: — ^• 



B. Moj 2 case 



The wavelength of the source field longer than Muj = 12 case and the wave effects become 
more significant. Fig. 15 shows spatial distribution of waves around the black hole. The 
source position is rs = 6M. 
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FIG. 15: Spatial distribution of $ on {z,x) plane for Mcj = 2. Left panel: Re[$]. Right panel: 
Im[$]. The point source is located at {z,x) = (— 6M, 0). 



The scattering amplitude also shows interference pattern. The distance between adjacent 
fringes becomes larger compared to Muj = 12 case. But the characteristic feature is the 
same. 
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FIG. 16: The scattering amplitude at robs = 20M. Mu = 2, rs = 6M. 

In this case, we can also observe a ring image for = corresponding to the unstable 
orbit. 
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FIG. 17: Images of black holes reconstructed from scattering waves. From the top left panel to the 
bottom right panel, the scattering angles are — 0, 7r/4, 7r/2, 37r/4, tt. 
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FIG. 18: The intensity distribution of images. Left panel: Oq = 0. Right panel: — ^^ 

Finally, we present the result for the source position rs = 2.5M. In this case, we can 
identify both the forward and the backward glories as ring images (Fig. 21). 
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FIG. 19: Spatial distribution of $ on {z,x) plane for Mcj = 2. Left panel: Re[$]. Right panel: 
Im[$]. The point source is located at {z,x) = (— 2.5M, 0). 
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FIG. 20: The scattering amplitude at robs = 20M. Mcj = 2,rs = 2.5M. 
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FIG. 21: Images of black holes reconstructed from scattering waves. From the top left to the 
bottom right panel, the scattering angles are — 0, 7r/4, 7r/2, 37r/4, tt. 
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FIG. 22: The intensity distribution of images. Left panel: — 0. Right panel: — ^^ 

V. SUMMARY AND DISCUSSION 

We investigated scattering of scalar waves by the Schwarzschild black hole. Our main 
purpose of this analysis is to obtain images of black holes using waves. For this purpose, 
we solved the scalar wave equation in the Schwarzschild spacetime numerically and recon- 
struct ecd images of black holes by the Fourier transformation of the scattering waves. For 
the forward and the backward direction of the scattering, we obtained ring images corre- 
sponding to the glory scattering by the black hole. In the geometric optics limit, these ring 
images are related to existence of the unstable orbit for null rays. As extension of analysis 
presented in this paper, wave scattering and image formation in the Kerr geometry is an 
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interesting subject to be tackled. For the Kerr black hole, due to dragging of the spacetime, 
incident waves can be amplified by the superradiance [TO]. By investigating images from 
scattered waves by the Kerr black hole, we expect to find out a new aspect and an interpre- 
tation of the superradiance in connection with the Penrose process. We will report on the 
analysis of the Kerr black hole case in our next publication. 

Another application of the analysis presented in this paper is related to observations of 
black hole shadows ^4[[T^. As the apparent angular sizes of black hole shadows are so small, 
the diffraction effect on images are crucial to resolve black hole shadows in observation using 
radio interferometer. For SgrA*, which is the black hole candidate at Galactic center, the 
apparent angular size of its shadow is estimated to be ~ 30/x arc seconds and this value is the 
largest among black hole candidates. For a sub-mm VLBI, a baseline length must be larger 
than 1000 km to resolve the shadow and this requirement shows the possibility to detect 
the black hole shadow of SgrA* using the present day technology of VLBI telescope. Thus, 
analysis of black hole shadows based on wave optics is an important thema for sucessful 
detection of shadows and detemination of black hole parameters via imaging of black holes. 
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